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ABSTRACT 

We  determine  the  masses  and  magnitude  difference  of  the  components  of  the  Hyades  spectroscopic  binary  92  Tauri. 
We  find  that  both  components  appear  to  be  less  massive  and/or  brighter  than  predicted  from  some  recent  evolutionary 
models.  The  rapid  rotation  and  unknown  rotational  inclination  of  both  components  introduce  uncertainty  in  their 
luminosities  and  colors,  but  not  enough  to  reconcile  both  of  them  with  the  evolutionary  models.  We  measured  the 
visual  orbit  with  the  Mark  III  optical  interferometer  and  the  Navy  Prototype  Optical  Interferometer  and  combined  it 
with  the  Hipparcos  proper-motion-based  parallax  to  find  a  total  system  mass  £A4  of  4.03  ±  0.20  A40.  We  also 
combined  our  visual  orbit  with  three  recent  spectroscopic  orbits  to  find  three  spectroscopically  based  estimates  of 
YiM  and  compared  these  to  the  EA4  from  the  visual  orbit  and  parallax.  We  chose  the  spectroscopic  orbit  that  agreed 
best  and  used  its  mass  ratio  to  estimate  individual  masses  AIa.b  of  2.15  ±  0.12  and  1.87  ±  0.11  AT-,.  From  the 
interferometry,  we  determine  Am  =  1.13  ±  0.05  mag  across  the  450-850  nm  band.  The  parallax  then  implies  ab¬ 
solute  V  magnitudes  M \  B  of  0.48  ±  0.05  and  1 .61  ±  0.06  mag.  If  the  components  are  rotating  near  breakup  veloc¬ 
ity  and  seen  nearly  pole-on,  the  true  luminosities  may  be  as  faint  as  1 .03  and  2.13  mag;  even  in  that  case,  however,  the 
secondary  is  too  blue  by  ~0.07  mag  in  B  —  V. 

Key  words:  binaries:  spectroscopic  —  binaries:  visual  —  stars:  individual  ( 9 2  Tauri)  — 
techniques:  interferometric 
Online  material:  OIF1TS  file 


1.  INTRODUCTION 

The  binary  92  Tauri,  the  brightest  member  of  the  Hyades, 
excited  some  interest  in  the  1990s  as  the  19  mas  visual  orbit  was 
resolved  by  interferometry  (Armstrong  et  al.  1991;  Pan  et  al. 
1992;  Hummel  &  Armstrong  1992),  while  the  spectral  lines  of 
the  rapidly  rotating  secondary  were  disentangled  from  those  of 
the  primary  (also  a  rapid  rotator)  to  yield  a  double-lined  spec¬ 
troscopic  orbit  by  Peterson  (1991,  hereafter  P91),  Peterson  et  al. 
(1993,  hereafter  PSL93),  Tomkin  et  al.  (1995,  hereafter  TPM95), 
and  Torres  et  al.  (1997,  hereafter  TSL97).  In  the  early  1990s,  the 
lure  of  92  Tau  was  the  orbital  parallax — and  thus  the  distance 
to  the  Hyades — to  be  derived  from  the  visual  and  spectroscopic 
elements. 

By  the  late  1990s,  the  Hipparcos  trigonometric  parallax  was 
available  (Perryman  et  al.  1998),  with  precision  comparable  to 


or  better  than  that  of  the  orbital  parallax.  In  200 1 ,  de  Bruijne  et  al. 
published  an  even  more  precise  proper-motion-based  parallax 
of  92  Tau  in  the  process  of  evaluating  the  parallaxes  of  a  large 
number  of  Hyades  members  and  determining  the  distance  to  the 
cluster.  As  a  result,  92  Tau  is  now  interesting  not  so  much  for  its 
distance  as  for  the  constraints  that  the  masses  and  luminosities  of 
its  components  impose  on  the  evolutionary  state  and  elemental 
abundances  of  the  Hyades.  In  particular,  its  components,  along 
with  5 1  Tau  A,  are  the  only  Hyades  stars  near  the  turnoff  point 
with  measured  masses. 

The  problems  with  using  the  masses  of  92  Tau  currently  in  the 
literature  are  that  the  visual  orbit,  based  on  data  from  the  Mark  III 
optical  interferometer,  has  been  reported  only  in  preliminary 
versions,  while  there  are  not  one,  but  three,  recent  spectroscopic 
orbit  calculations  with  somewhat  discordant  results.  The  spectro¬ 
scopic  results  and  the  implied  masses,  distances,  and  magnitudes 
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TABLE  1 

Recent  Spectroscopic  Results 


Parameter  PSL93/P913  TPM95b  TSL97b 


Spectroscopy 

Ka  (km  s-1) .  31.2  ±  0.9  29.3  ±  0.3  33.18  ±  0.39 

Kb  (km  s-1) .  42.7  ±  2.4  38.0  ±  2.8  38  ±  2 

va  sin i  (km  s_1) .  65  70  ±  5  70 

vb  sin  i  ( km  s-1) .  170  90  ±  20  110  ±  4 

q .  0.73  ±  0.05  0.77  ±  0.06  0.873  ±  0.048 


With  Visual  Elements  Added 


^orb  (mas) .  21.86  ±  0.90  22.7  ±  1.1  21.22  ±  0.76 

D  (pc) .  45.7  ±  1.9  44.1  ±  2.2  47.10  ±  1.7 

m  -  M  (mag) .  3.302  ±  0.089  . . .  3.37  ±  0.08 

Ma(M0) .  2.48  ±  0.37  2.1  ±  0.3  2.42  ±  0.30 

Mb(M0) .  1.80  ±0.17  1.6  ±0.2  2.11  ±  0.17 


a  As  calculated  from  PSL93  using  correct  AvR  data  from  P9 1  (see  text)  and 
visual  orbital  elements  from  this  work. 

b  As  published  by  TPM95  and  TSL97  using  visual  orbital  elements  from 
Pan  et  al.  (1992). 

are  summarized  in  Table  1  (with  corrections  to  the  PSL93  re¬ 
sults;  see  §  2.1).  The  de  Bruijne  et  al.  (2001)  parallax  and  the 
visual  interferometric  orbit  we  present  here  significantly  clarify 
the  situation,  as  we  show  below. 

Our  visual  orbit  is  based  on  30  nights  of  data  from  the  Mark  III 
optical  interferometer,1  which  ceased  operations  at  the  end  of 
1992,  and  four  nights  of  data  from  the  Navy  Prototype  Optical 
Interferometer2  (NPOI),  which  has  operated  at  the  Lowell  Ob¬ 
servatory  site  on  Anderson  Mesa,  Arizona,  since  1 995.  The  NPOI 
data  improve  the  sampling  near  periastron  and  include  closure 
phases,  which  permit  us  to  resolve  the  1 80°  ambiguity  in  the  po¬ 
sition  angle  of  the  line  of  nodes.  The  resulting  orbit  is  consistent 
with,  but  considerably  more  precise  than,  the  preliminary  or¬ 
bits  published  by  Armstrong  et  al.  ( 1 99 1 ),  Pan  et  al.  ( 1 992),  and 
Hummel  &  Armstrong  (1992)  based  on  subsets  of  the  Mark  III 
data. 

Besides  presenting  the  visual  orbit,  the  purpose  of  this  paper  is 
to  infer  the  best  component-mass  estimates  by  combining  the  best 
spectroscopic  orbit  with  our  visual  orbit.  To  select  the  best  spec¬ 
troscopic  orbit,  we  first  use  the  fact  that  the  total  system  mass  can 
be  determined  by  combining  the  visual  orbit  with  either  a  spec¬ 
troscopic  orbit  or  the  parallax.  We  then  choose  the  spectroscopic 
orbit  that  produces  a  total-mass  estimate  that  agrees  best  with  the 
parallax-based  total  mass.  Finally,  we  calculate  component  masses 
from  the  ratio  of  the  velocity  amplitudes.  Because  the  parallax 
constrains  the  total  mass,  the  uncertainties  of  the  component 
masses  are  considerably  smaller  than  those  obtained  from  the 
combination  of  the  visual  and  spectroscopic  orbits  alone.  The 
interferometry  also  produces  a  measurement  of  magnitude  dif¬ 
ferences  across  the  observed  wavelength  range. 

As  we  show  below,  the  component  masses,  luminosities,  and 
colors  present  some  problems  for  evolutionary  models  that  have 
been  applied  to  the  Hyades.  Determining  the  luminosities  and 
the  colors  of  the  components  is  complicated  by  their  rapid  ro¬ 


1  The  Mark  III  optical  interferometer,  located  on  Mount  Wilson,  was  jointly 
operated  by  the  Naval  Research  Laboratory  and  the  US  Naval  Observatory,  with 
funding  from  the  Office  of  Naval  Research  and  from  the  Oceanographer  of  the 
Navy. 

2  The  Navy  Prototype  Optical  Interferometer  is  a  joint  project  of  the  Naval 

Research  Laboratory  and  the  US  Naval  Observatory,  in  cooperation  with  Lowell 
Observatory. 


tation.  The  recent  discussion  by  Peterson  et  al.  (2006)  of  the  ef¬ 
fects  of  rotation  on  Altair  demonstrates  that  the  observed  colors 
and  fluxes  of  a  rapid  rotator  can  differ  significantly  from  those  of 
an  equivalent  nonrotating  star,  depending  on  the  inclination  of 
the  rotation  axis  to  the  line  of  sight.  We  have  calculated  correc¬ 
tions  to  B  —  V  and  My  based  on  the  range  of  possible  inclina¬ 
tions.  Even  with  the  range  of  possible  corrections,  problems  with 
the  evolutionary  models  remain. 

We  summarize  the  observational  background  for  02  Tau  in  §  2. 
We  then  describe  the  observations  and  analysis  on  which  our 
visual  orbit  is  based  (§§  3  and  4),  present  the  resulting  masses  in 
§  5,  and  compare  them  with  evolutionary  models  in  §  6. 

2.  BACKGROUND 
2.1.  Spectroscopic  and  Visual  Orbits 

The  binary  nature  of  02  Tau  (=  HR  1412  =  HD  283 19  =  vB  72  = 
HIP  20894;  a,  6  [J2000.0]  =  4h28m39?3, 15°52'16"4;  mv  =  3.41; 
B  —  V  =  0.18;  Sp  =  A7  III  [it  is  a  standard  for  this  spectroscopic 
type,  despite  its  duplicity])  was  first  noted  by  Moore  (1908)  and 
Frost  (1909).  Plaskett  (1915),  Petrie  (1940),  and  Ebbighausen 
(1959)  measured  the  spectroscopic  orbit  of  the  primary  but  failed 
to  detect  the  secondary,  despite  the  fact  that  the  magnitude  dif¬ 
ference  is  only  slightly  greater  than  1  mag. 

The  secondary  was  first  observed  in  lunar  occultation  obser¬ 
vations  (Evans  &  Edwards  1980;  Peterson  et  al.  1981),  and  the 
difference  A vr  between  the  radial  velocities  of  the  components 
was  measured  at  periastron  using  cross-correlation  techniques 
by  P9 1 .  The  reason  for  the  earlier  lack  of  detection  is  the  large 
rotational  velocity  of  the  secondary  («>b  sin  i  is  given  variously 
between  90  and  170  km  s_1);  its  lines  are  broad  and  diffuse,  and 
always  blended  with  those  of  the  primary. 

The  Mark  III  optical  interferometer  resolved  the  visual  orbit 
in  observations  during  1989-1991,  with  preliminary  results  re¬ 
ported  by  Armstrong  etal.  (1991),  Pan  etal.  (1992),  and  Hummel 
&  Armstrong  (1992).  During  the  1990s,  three  groups — PSL93, 
followed  at  2  year  intervals  by  TPM95  and  TSL97 — obtained 
spectroscopic  data  and  combined  the  velocity  amplitudes  with 
the  Pan  et  al.  ( 1 992)  results  to  derive  masses  and  the  distance. 
However,  the  velocity  amplitudes  from  the  three  groups  are  some¬ 
what  discordant,  apparently  due  to  their  differing  approaches  to 
the  line-blending  problem.  Their  results — velocity  amplitudes 
K,  rotational  velocities  v  sin  i,  and  mass  ratios  q — are  summa¬ 
rized  in  the  first  half  of  Table  1 . 

PSL93  used  synthetic  spectra  based  on  Kurucz  (1991)  model 
atmospheres  to  separate  the  contributions  of  the  two  compo¬ 
nents.  They  found  a  maximum  cross-correlation  between  syn¬ 
thetic  and  observed  spectra  for  «>a,b  sin  i  =  65  and  170  km  s  ’, 
but  the  radial  velocity  of  the  secondary  depended  so  sensitively 
on  its  rotational  velocity  that  they  could  not  derive  a  numerically 
stable  value  for  its  amplitude.  Nevertheless,  they  were  able  to 
remove  the  effects  of  the  secondary  from  the  velocities  of  the 
primary.  They  then  used  velocity  differences  from  P91  with  the 
improved  velocities  of  the  primary  to  produce  a  spectroscopic 
orbit  that  included  both  components. 

Unfortunately,  the  KB  results  reported  by  PSL93  were  in  er¬ 
ror  because  of  a  mistake  in  transcribing  the  Avr  data  from  P91 
(D.  M.  Peterson  1994,  private  communication).  We  have  recalcu¬ 
lated  Ka  b  and  q  using  velocities  of  the  primary  given  by  PSL93 
and  the  Avr  data  of  P9 1 ;  these  are  the  results  that  are  shown  in 
Table  1. 

TPM95  cross-correlated  their  6 2  Tau  spectra  with  the  spectra 
of  radial  velocity  standards  to  get  the  velocity  curve  of  the  pri¬ 
mary.  They  then  subtracted  the  velocity-shifted  spectrum  of  HR 
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1427,  an  A6  IV  Hyades  member  with  a  projected  rotational 
velocity  approximately  equal  to  that  of  62  Tau  A,  from  their  data 
to  produce  spectra  of  92  Tau  B  alone.  From  these  spectra  they 
obtained  the  velocity  curve  of  the  secondary  and  estimated  that 
vb  sin  i  =  90  km  s_1  by  comparison  with  the  broadened  spec¬ 
trum  of  a  sharp-lined  star. 

TSL97,  like  PSL93,  compared  their  spectra  to  synthetic  spec¬ 
tra.  The  synthetic  spectra  were  composed  of  templates  with  T^f  = 
8250  K  for  both  stars,  log  g  =  4.0  and  4.5  for  the  primary  and 
secondary,  respectively,  and  va  sin  i  =  70  km  s_1,  with  t>B  sin  i 
varying  between  90  and  170  km  s_1.  Because  they  were  unable 
to  obtain  reliable  secondary  radial  velocities  from  individual  spec¬ 
tra,  they  took  a  global  approach  and  found  that  Kb  =  38.0kms_1 
and  vb  sin  /  =  110  km  s~'  gave  the  best  average  correlation  over 
all  exposures. 

The  orbital  parallaxes  zz7orb,  distances  D,  and  masses  M  de¬ 
rived  by  these  three  groups  are  summarized  in  the  second  half  of 
Table  1  (with  corrected  results  from  PSL93).  The  distances  are 
consistent  within  the  uncertainties,  as  are  the  masses  of  the  pri¬ 
maries.  The  secondary  masses  are  slightly  more  discordant;  the 
TPM95  and  TSL97  ranges  do  not  overlap,  although  the  discrep¬ 
ancy  is  not  large. 

These  distance  and  mass  estimates  have  been  cited  by  several 
authors.  Narayan  &  Gould  (1999)  used  92  Tau  with  two  other  bi¬ 
naries  to  check  for  systematic  errors  in  the  Hipparcos  parallaxes. 
Lastennet  et  al.  (1999)  used  92  Tau,  also  with  two  other  binaries, 
to  evaluate  stellar  evolution  codes.  Lebreton  et  al.  (2001)  used 
92  Tau  with  four  other  binaries  to  construct  an  empirical  mass- 
luminosity  function  and  to  assess  the  age  and  helium  content  of 
the  cluster. 

2.2.  Astrometry 

In  1997  the  Hipparcos  trigonometric  parallaxes  became  avail¬ 
able.  For#2  Tau,  w  =  21.89  ±  0.83  mas  (Perryman  et  al.  1997), 
as  precise  as  the  orbital  parallaxes.  Because  the  center  of  light 
moves  over  a  range  projected  on  the  sky  of  «4.5  mas,  the  fact 
that  the  Hipparcos  reduction  treated  92  Tau  as  a  single  star  is  a 
possible  reason  for  concern.  It  is  likely,  however,  that  the  only  ef¬ 
fect  is  reduced  precision,  for  several  reasons.  First,  the  residuals 
of  the  Hipparcos  solution  binned  with  the  orbital  period  show 
only  a  very  weak  correlation  with  orbital  phase,  if  any.  Second, 
the  orbital  period  is  short  compared  to  the  length  of  the  Hipparcos 
mission.  Finally,  the  large  eccentricity  results  in  little  motion  of 
the  center  of  light  during  most  of  the  orbital  period. 

Dravins  et  al.  (1997)  used  Hipparcos  parallaxes  and  proper 
motions  to  determine  the  distance  to  the  Hyades  and  parallaxes 
for  its  members,  although  they  did  not  list  individual  parallaxes. 
They  noted  that  the  proper-motion-based  parallaxes  were  about 
twice  as  precise  as  the  trigonometric  parallaxes.  De  Bruijne  et  al. 
(2001)  used  both  Hipparcos  and  Tycho  data  to  determine  proper- 
motion  parallaxes.  For  92  Tau,  they  give  a  Hipparcos  parallax  of 
22.24  ±  0.36  mas  and  a  Tycho  parallax  of  22.35  ±  0.36  mas. 

2.3.  b  Scuti  Pulsations 

In  a  survey  of  Hyades  A-F  stars,  Horan  (1977,  1979)  discov¬ 
ered  that  92  Tau  is  a  b  Scuti  variable.  The  pulsations  were  con¬ 
firmed,  and  more  frequencies  measured,  by  Breger  et  al.  (1987, 
1989);  in  addition,  Breger  et  al.  (1989)  determined  that  it  is  the 
primary  that  is  pulsating.  New  photometry  by  Breger  et  al.  (2002) 
may  show  pulsations  of  the  secondary  as  well.  From  the  mass 
estimate  published  by  TSL97,  evolutionary  modeling,  and  the 
observed  pulsation  frequencies,  Breger  et  al.  (2002)  determined 
that  92  Tau  A  is  close  to  the  blue  edge  of  the  instability  strip  and 
that  Teff  Si  7800  K.  The  lower  mass  estimate  in  our  work  moves 


92  Tau  A  away  from  the  edge  of  the  strip,  relaxing  the  constraint 
on  reff. 

3.  OBSERVATIONS  AND  CALIBRATION 

We  observed  92  Tau  during  1989-1991  with  the  Mark  III 
optical  interferometer  and  obtained  30  nights  of  useful  data,  with 
between  three  and  42  scans  per  night  in  three  spectral  channels 
at  800,  550,  and  450  nm  in  1989,  and  at  800,  550,  and  500  nm 
thereafter.  Filter  widths  were  22  nm  for  800  nm  and  25  nm  for  the 
other  wavelengths.  The  baseline  lengths  used  varied  between  8.2 
and  3 1 .5  m.  A  Mark  III  scan  consisted  of  one  75  s  stellar  obser¬ 
vation  in  three  bands  plus  a  5  s  background/dark  count  mea¬ 
surement.  Descriptions  of  the  Mark  III  in  some  detail  are  given 
by  Shao  et  al.  (1988)  and  Mozurkewich  et  al.  (1991),  and  more 
briefly  in  the  context  of  binary  observations  by,  e.g.,  Armstrong 
etal.  (1992). 

We  also  observed  92  Tau  in  1997  and  1998  with  the  NPOl  and 
obtained  four  nights  of  useful  data,  with  between  four  and  seven 
scans  per  night  in  19  spectral  channels  covering  850-520  nm. 
We  used  a  configuration  consisting  of  the  center,  east,  and  west 
astrometric  stations,  with  baseline  lengths  of  37.5,  22.2,  and 
18.9  m.  These  scans  consisted  of  one  90  s  observation  plus  a  5  s 
background/dark  count  measurement.  A  description  ofthe  NPOl 
is  given  by  Armstrong  et  al.  (1998);  binary  observations  with  the 
NPOl  are  described  by,  e.g.,  Benson  et  al.  (1997)  and  Hummel 
et  al.  (1998).  The  calibration  of  these  data  followed  the  methods 
we  have  used  in  previous  work  with  Mark  III  and  NPOl  data, 
e.g.,  Mozurkewich  et  al.  (1991)  and  Armstrong  et  al.  (1992)  for 
the  Mark  III,  and  Benson  et  al.  (1997)  and  Hummel  et  al.  (1998) 
for  the  NPOl. 

The  contributions  of  the  two  instruments  to  these  results  are 
about  equal:  1052  visibilities  in  30  nights  from  the  Mark  III,  plus 
785  visibilities  and  122  closure  phases  in  four  nights  from  the 
NPOl.  The  separations  from  the  four  NPOl  nights  are  in  general 
more  precise  because  of  the  greater  number  of  visibilities  per 
night.  Table  2  lists  the  observations.  The  92  Tau  calibrated  vis¬ 
ibilities  and  closure  phases,  along  with  their  uncertainties,  are 
available  in  the  electronic  edition  of  the  Journal  in  the  FITS- 
based  Optical  Interferometry  Data  Exchange  format  (OIFITS) 
described  by  Pauls  et  al.  (2005). 3 

4.  DATA  REDUCTION  AND  ANALYSIS 

The  data  analysis  method  used  here,  described  in  detail  by 
Hummel  et  al.  (1993,  1998),  consists  of  using  the  Marquardt- 
Levenberg  method  (e.g.,  Press  et  al.  1992)  to  fit  a  set  of  com¬ 
ponent  parameters  and  orbital  parameters  directly  to  the  entire 
set  of  measured  visibilities  V2(A,  t ).  The  first  92  Tau  results  using 
this  technique  were  reported  by  Hummel  &  Armstrong  (1992). 
Fitting  the  parameters  directly  to  the  V2  data,  bypassing  the  stage 
of  fitting  (y o,9)  values  to  each  night’s  visibilities  used  by  Armstrong 
et  al.  (1991)  and  Pan  et  al.  ( 1992),  allows  us  to  use  nights  on  which 
we  took  only  a  few  scans  and  nights  close  to  periastron  during 
which  the  change  in  9  was  not  negligible.  Fitting  to  the  entire  V2 
data  set  also  gives  us  a  more  direct  estimate  of  the  parameter 
uncertainties. 

The  component  parameters  in  the  fit  include  the  disk  diame¬ 
ters  (7a,  b  and  the  magnitude  difference  parameterized  by  the  dif¬ 
ferences  at  500, 550,  700,  and  800  nm.  (Magnitude  differences  at 
other  wavelengths  were  estimated  by  a  spline  interpolation. )  The 
orbital  parameters  are  the  standard  seven:  the  angular  semimajor 
axis  of  the  relative  orbit  a",  eccentricity  e,  period  P.  epoch  of 


3  Defined  at  http://www.mrao.cam.ac.uk/~jsyl001/exchange. 
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TABLE  2 

Observation  and  Result  Log  for  02  Tau 

HJD  BL  p  9  \0  -  C\  a0-c 

(+2,440,000)  Orbital  Phase  (m)  Number  of  Scans  (mas)  (deg)  (mas)  (mas) 

(1)  (2)  (3)  (4)  (5)  (6)  (7)  (8) 

1989  (Mark  III) 

7772.967 .  0.195  31.5  6  21.24  7.88  0.17  0.30 

7782.952  .  0.266  8.2  3  24.38  15.16  0.44  1.92 

7790.961 .  0.323  23.1  4  26.02  23.62  1.02  0.97 

7791.954 .  0.330  23.1  8  25.94  22.37  0.16  0.56 

7804.952  .  0.422  31.5  8  26.34  28.88  0.51  0.54 

7805.952  .  0.429  31.5  7  26.67  30.15  0.06  0.46 

7806.947 .  0.436  31.5  7  27.44  33.32  1.39  1.10 

7813.928 .  0.486  8.2  9  26.38  34.93  0.01  0.81 

7818.933 .  0.521  15.2  6  26.18  38.83  0.46  0.43 

7826.945 .  0.578  15.2  10  25.11  43.33  0.25  0.35 

7827.966  .  0.585  15.2  6  25.45  44.61  0.82  1.26 

7831.907 .  0.613  15.2  12  24.03  46.47  0.004  0.34 

7832.911  .  0.621  27.0  20  24.03  47.62  0.26  0.13 

7833.822 .  0.627  27.0  3  22.83  44.50  1.61  1.48 

7836.827  .  0.648  23.3  13  23.22  50.74  0.30  0.32 

7865.805 .  0.854  8.2  9  13.92  84.87  0.60  0.42 

1990  (Mark  III) 

8148.966  .  2.866  15.4  8  14.14  88.52  0.23  0.12 

8150.917 .  2.880  27.0  3  14.00  93.85  0.82  0.67 

8187.803  .  3.142  23.6  7  17.77  359.06  0.45  0.57 

8189.953  .  3.158  27.6  38  18.81  2.80  0.07  0.040 

8192.922 .  3.179  31.5  42  20.23  6.12  0.03  0.11 

8233.833 .  3.470  31.5  17  26.53  33.83  0.08  0.08 

8234.827  .  3.477  31.5  20  26.49  34.39  0.07  0.17 

1991  (Mark  III) 

8517.951 .  5.488  27.0  13  26.26  35.05  0.10  0.12 

8548.927  .  5.709  15.2  18  20.86  57.32  0.14  0.14 

8549.928  .  5.716  15.2  17  20.56  58.20  0.16  0.28 

8562.911  .  5.808  15.2  15  17.08  73.88  0.30  0.23 

8593.872 .  6.028  27.6  15  4.62  250.21  1.16  0.30 

8594.793 .  6.034  27.6  12  4.69  274.87  0.88  0.22 

8597.728 .  6.055  15.2  11  6.89  323.65  0.24  0.12 

1997  (NPOI) 

10459.722 .  19.286  37.5,  22.2,  18.9  6  24.93  18.29  0.06  0.067 

10736.803  .  21.255  37.5,  22.2,  18.9  7  24.15  15.17  0.17  0.15 

1998  (NPOI) 


10836.680 .  21.965  37.5,  22.2,  18.9  7  8.92  134.44  0.14  0.034 

11131.849 .  24.063  37.5,22.2,18.9  4  9.16  333.92  1.10  0.33 


Notes. — Col.  (1):  Heliocentric  Julian  Date  of  observation.  Col.  (2):  Orbital  phase  relative  to  the  epoch  determined  in  this  work  and  the  period  from  TPM95. 
Col.  (3):  Baseline  length(s).  Col.  (4):  Number  of  scans  available  for  that  night.  The  orbital  solutions  used  the  fringe  visibilities  directly,  rather  than  the  intermediate 
solutions  for  each  night’s  separation  p  and  position  angle  6\  however,  to  indicate  the  quality  of  the  data,  the  remaining  columns  give  the  results  for  fits  to  individual 
nights.  Cols.  (5)  and  (6):  Best-fit  separation  p  and  position  angle  6  for  each  night.  Col.  (7):  Distance  from  the  best-fit  position  to  the  position  predicted  by  the  best-fit 
orbit.  Col.  (8):  Distance  from  the  best-fit  position  to  the  edge  of  the  68%  confidence  region,  as  measured  along  the  O  —  C  vector. 

periastron  T,  orbital  inclination  iorh ,  argument  of  periastron  u>,  (Peterson  et  al.  1981;  PSL93).  We  then  applied  the  formulas 

and  position  angle  of  the  line  of  nodes  Cl.  given  in  Mozurkewich  et  al.  (1991)  and  obtained  uniform-disk 

In  the  case  of  61  Tau,  we  did  not  fit  for  the  component  di-  diameter  estimates  at  800  nm  of  0.76  and  0.42  mas.  Because  both 

ameters  dA  B .  Instead,  we  estimated  the  diameters  from  the  in-  angular  diameters  are  below  the  resolutions  of  the  Mark  III  and 

dividual  R  —  I  colors  and  apparent  V  magnitudes.  We  inferred  the  NPOI  configuration  used  (the  smallest  fringe  spacing,  at 

the  apparent  magnitudes  of  the  components  from  the  system  520  nm  on  the  37.5  m  NPOI  baseline,  was  ss2.8  mas),  errors  in 

magnitude  and  the  magnitude  difference  Am  from  our  data.  As  a  the  adopted  diameters  have  a  negligible  effect  on  the  estimates 

starting  point,  we  took  the  individual  R  —  I  colors  to  be  equal  to  for  other  system  parameters. 

the  combined  color  (0.09  mag;  Johnson  et  al.  1966)  based  on  the  We  adopted/3  =  140.72816  ±  0.00093  days  from  TSL97,  since 

observation  that  the  colors  of  the  components  are  nearly  equal  it  is  based  on  observations  spanning  90  years  rather  than  the  9  years 
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covered  by  our  observations.  Then,  with  P  and  dA  B  selected,  we 
determined  the  other  six  orbital  elements,  redetermined  Am  at 
four  wavelengths,  and  repeated  the  orbit  fitting  after  a  slight  ad¬ 
justment  of  dA  B.  This  procedure  converged  after  one  iteration. 

The  value  of  xl  at  the  best  solution  was  2.45.  The  most  discre¬ 
pant  point  in  the  orbit,  in  units  of  ao-c,  is  at  JD  2,448,594.793 
[(p,  9)  =  (4.69  mas,  274?87)];  when  it  is  removed,  drops 
to  2.39  and  the  best  solution  shifts  slightly.  In  the  absence  of  a 
clear  reason  to  exclude  it,  we  have  kept  that  night  in  the  solution. 
(The  point  at  JD  2,450,836.680  [(p,  9)  =  (8.92  mas,  134?44)] 
is  more  discrepant  until  the  wavelength  uncertainty  of  the  NPOl, 
discussed  below,  is  taken  into  account.) 

We  adopted  the  conservative  assumption  that  differs  from 
unity  due  to  underestimation  of  the  observational  uncertainties. 
Using  increased  observational  uncertainties,  we  estimated  the  er¬ 
rors  in  the  orbital  elements  from  the  curvature  of  the  \2  hyper¬ 
surface.  We  then  checked  these  error  estimates  by  confirming  that 
perturbing  the  solution  by  1  a  in  one  parameter  and  then  fitting 
the  remaining  parameters  produced  an  increase  in  of  one  part 
in  1949,  the  number  of  degrees  of  freedom  in  the  fit. 

Our  final  step  was  to  produce  (p,  9)  values  and  68%  confi¬ 
dence  error  ellipses  for  each  night  in  order  to  illustrate  the  orbit. 
The  distances  from  the  nightly  positions  predicted  by  the  best-fit 
orbit  and  those  from  the  single-night  fits  are  shown  in  column  (7) 
of  Table  2.  Column  (8)  lists  the  distance  from  each  single-night 
fit  to  the  edge  of  the  68%  confidence  region,  as  measured  along 
the  line  toward  the  predicted  position  for  that  night. 

Included  in  the  error  estimate  for  the  NPOl  nights  is  a  poten¬ 
tial  error  in  the  angular  scale  due  to  the  uncertainty  in  the  wave¬ 
length  scale.  In  the  NPOl  spectrometers,  the  dispersed  light  is 
focused  on  an  array  of  lenslets  that  define  the  spectral  channels. 
The  lenslet  arrays  are  aligned  using  the  HeNe  laser  that  is  also 
used  for  aligning  the  beam-combiner  table,  delay  lines,  and  feed 
system.  The  alignment  is  stable  and  accurate  at  the  :S0.5%  level, 
which  we  adopt  as  our  wavelength  uncertainty  for  the  NPOl 
data.  Note  that  the  Mark  III  used  filters  with  well-measured  and 
stable  passbands  and  that  the  shift  in  2eff  f°r  the  Mark  III  due  to 
the  spectral  shape  of  an  A  star  is  less  than  one  part  in  3500. 

Because  the  angular  scale  is  proportional  to  1/2,  this  uncertainty 
propagates  directly  into  an  uncertainty  in  the  binary  separation  and 
in  a".  The  formal  precision  of  a"  is  very  high:  18.796  ±  0.011  mas, 
or  one  part  in  1700.  Because  about  half  of  these  data  were  taken 
with  the  Mark  III,  for  which  the  wavelengths  were  well  measured 
and  stable,  we  adopt  0.3%  as  the  uncertainty  due  to  the  wave¬ 
length  scale.  This  value  dominates  the  error  budget  and  results  in 
an  uncertainty  in  a"  of  0.056  mas. 

5.  RESULTS  AND  DISCUSSION 

The  results  of  the  visual  orbit  fitting  are  shown  in  Figure  1  and 
listed  in  Table  3,  along  with  the  elements  from  Armstrong  et  al. 
(1991),  Pan  et  al.  (1992),  and  Hummel  &  Armstrong  (1992).  Our 
orbit  is  consistent  with  the  earlier,  preliminary  elements  but  has 
considerably  higher  precision,  which  we  attribute  to  the  method 
of  directly  fitting  the  V2  data,  the  larger  data  set  used,  and  im¬ 
proved  coverage  near  periastron. 

With  the  visual  orbit  in  hand,  determining  the  component 
masses  proceeds  in  four  steps.  First,  we  determine  the  total  mass 
of  the  system  from  a",  w,  and  P,  a  result  that  is  independent  of 
Ka  B .  Second,  we  determine  the  total  mass  implied  by  e,  ;or b ,  and 
the  Ka  b  values  of  PSL93  (as  corrected  from  P91),  TPM95,  or 
TSL97.  These  masses  are  independent  of  w.  Third,  we  choose 
the  AZAb  set  that  gives  the  best  agreement  with  the  first  step. 
Finally,  we  divide  the  total  mass  between  the  components  using 
q,  the  ratio  of  KB  to  KA. 


A-B 


Fig.  1 . — Visual  orbit  of  92  Tau.  The  ellipse  shows  the  best-fit  orbit  as  de- 
termined  by  a  direct  fit  to  the  interferometric  visibilities  (see  text).  The  center 
of  each  error  ellipse  marks  the  best-fit  separation  for  one  night’s  observations 
and  is  connected  to  the  corresponding  position  on  the  best-fit  orbit,  while  its 
size  shows  the  67%  confidence  region  for  that  night’s  data.  Some  error  ellipses 
are  nearly  too  small  to  see  in  this  figure,  e.g.,  the  point  near  (A a,  A<5)  =  (6.4, 
—6.2).  The  line  from  (0,  0)  shows  the  position  angle  of  periastron. 


The  total  mass  of  the  system,  E M ,  can  be  calculated  directly 
from  the  parallax,  the  semimajor  axis,  and  the  period  inserted 
into  Kepler’s  third  law: 


E  M  +  Mr  = 


(a" /mf 

—p2 

yr 


M 


©’ 


(1) 


where  the  subscript  w  indicates  that  the  total  mass  is  derived  from 
the  parallax  and  PyI  is  of  course  the  period  in  years.  From  Table  3 
of  de  Bruijne  et  al.  (2001),  we  adopt  zu  =  22.30  ±0.36  mas,  the 
average  of  their  parallaxes  from  the  Hipparcos  and  the  Tycho 
proper  motions,  which  agree  very  well  with  one  another.  From 
our  visual  orbit,  we  adopt  a"  =  18.796  ±  0.056  mas,  and  from 
TSL97  we  adopt  Pyr  =  0.38529  yr  with  negligible  uncertainty. 
The  result  is  E_MOT  =  4.03  ±  0.20  AJ  0,  which  appears  in  the 
first  line  of  Table  4.  The  precision  of  this  result,  «6%,  is  domi¬ 
nated  by  the  precision  of  zu. 

As  we  discuss  below,  this  E M  w  is  smaller  than  the  expected 
value,  given  the  age  and  metallicity  of  the  Hyades  and  the  par¬ 
allax  and  apparent  magnitudes  of  92  Tau.  Note  that  using  the  Pan 
etal.  (1992)  value  of  18.6  ±  0.2  mas  for  a"  gives  an  even  smaller 
value, 3.91  ±  0.25  A40,  for  EA4ro.  On  the  other  hand,  using  the 
Perryman  et  al.  (1997)  Hipparcos  trigonometric  parallax  with 
our  a"  produces  a  larger  value,  4.26  ±  0.49  M  0 ,  with  a  consid¬ 
erably  larger  uncertainty. 

We  can  also  calculate  the  total  mass  from  the  velocity  ampli¬ 
tudes  Ka  b,  eccentricity  e,  and  orbital  inclination  forb: 


A/ AT  ^  —  AT  a  H-  AT  g 


1.036x  10  1Pi  ,TJr  „  ,3/i  2\3/2 


sin3  /0rb 


'-(KA+KB)3(l~e2)yi  (2) 
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TABLE  3 

Visual  Orbital  Elements  for  62  Tau 


Element 

A91a 

Pan92b 

H92c 

This  Work 

a  (mas) . 

18.6 

18.6  ±  0.2 

18.7 

18.796  ±  0.056 

e . 

0.73 

0.735 

0.73725  ±  0.00036 

P  (days) . 

140.75 

140.72816  (TSL97) 

T  -  2,440,000  (JD) ... 

7748.2 

7748.394 

7748.576  ±  0.012 

'orb  (deg) . 

46.7 

46.2  ±  1.0 

47.3 

47.61  ±  0.09 

w  (deg) . 

54 

57.7  ±  1.7 

54.8 

55.40  ±  0.06 

B  (deg) . 

174.3 

171.2  ±  1.8 

174.2 

173.73  ±  0.07 

Am450mn  (mag) . 

1.16  ±  0.13 

1.0 

A  (mag) . 

1.05  ±  0.05 

Az«550nm  (mag) . 

0.5 

1.18  ±  0.10 

1.0 

1.12  ±  0.03 

Amyoonm  (mag) . 

1.16  ±  0.03 

Arn800nm  (mag) . 

1.15  ±  0.07 

1.2 

1.14  ±  0.05 

a  Armstrong  et  al.  (1991). 
b  Pan  et  al.  (1992). 
c  Hummel  &  Armstrong  (1992). 


(e.g.,  Heintz  1978,  §  32),  where  the  subscript  K  indicates  that 
the  sum  is  derived  from  the  velocity  amplitudes,  Pj  is  the  pe¬ 
riod  in  days,  and  KA  B  are  in  kilometers  per  second. 

From  the  three  sets  of  KA  B  and  their  uncertainties,  combined 
with  our  values  for  ;orb  and  e  and  their  uncertainties,  plus  the 
period  of  TSL97  once  again,  we  calculated  the  three  values  for 
EA4  k  that  are  given  in  the  last  three  lines  of  Table  4.  These 
EA4  k  values  differ  from  the  sums  of  the  masses  in  Table  1  from 
previous  work  because  of  the  updated  values  for  zorb  and  e.  The 
E M  k  are  less  precise  than  E.M  w  because  of  the  difficulty  in 
measuring  the  velocities  in  this  system. 

The  best  agreement  between  E M  k  and  E M  w  comes  from  us¬ 
ing  the  TSL97  values  for KA  B ;  values  fromPSL93/P91  and  TPM95 
give  somewhat  discrepant  results.  Therefore,  we  use  q,  the  ratio 
of  the  velocity  amplitudes  from  TSL97,  to  allocate  EA4  between 
the  components.  We  find  {MA,  MB)  =  (2.15,  1.87)  A40. 

The  component  masses  are  correlated,  since  they  are  con¬ 
strained  by  EA4  and  q.  We  calculated  the  joint  probability 
V{MA,  M  b)  as  a  function  of  the  two  masses.  In  Figure  2  {left), 
we  show  the  boundaries  of  the  50%,  68%,  and  90%  confidence 
regions  in  this  joint  {MA,  M  b)  space.  In  particular,  note  that  the 
highest  total  mass  in  the  68%  confidence  region  is  at  M  A  = 
2.31  M  .  ,V( H  2.01  A/.. 

Our  interferometric  solution  also  includes  the  magnitude 
differences  A  m{X),  which  are  consistent  with,  but  more  precise 
than,  the  earlier  results.  In  principle,  calculating  the  components’ 
absolute  V  magnitudes  from  the  total  mv,  A  my,  and  parallax  is 
straightforward, yielding {MyAl  My,B)  =  { 0.48,  1.61)mag.  The 
correlation  between  these  magnitudes  is  shown  in  Figure  2  {right). 
The  lowest  luminosity  extreme  of  the  68%  confidence  region  lies 


TABLE  4 

Total  Mass  of  82  Tauri 


£  M 

Parameter  (A4  Q) 


(from  visual  orbit  plus  astrometry): 

a"  plus  zu  from  de  Bruijne  et  al.  (2001) .  4.03  ±  0.20 

TiM.k  (from  visual  and  spectroscopic  orbits): 

e  and  zorb  plus  KA$  from  P91/PSL93 .  4.50  ±  0.47 

e  and  zorb  plus  KAB  from  TPM95 .  3.40  ±  0.43 

e  and  zorb  plus  KAq  from  TSL97 .  4.02  ±  0.35 


at  My,  a  =  0-54  mag  and  My, b  =  1-68  mag.  In  practice,  as  we 
discuss  toward  the  end  of  the  next  section,  the  rotation  rate  and 
inclination  affect  the  apparent  magnitude  and  color  of  each  com¬ 
ponent,  which  complicates  the  comparison  with  evolutionary 
models.  The  component  masses  and  magnitudes  are  summa¬ 
rized  in  Table  5. 

6.  COMPARISON  WITH  EVOLUTIONARY  MODELS 
6.1.  Observed  Magnitudes  and  Colors 

The  masses  and  absolute  magnitudes  we  derive  here  pose  a 
problem  in  the  context  of  some  recent  sets  of  isochrones.  If  we 
accept  the  current  estimates  of  the  helium  content  (7  ss  0.28) 
and metallicity  (the  spectroscopic  [Fe/H]  is«0.14,  orZ  «  0.026) 
of  the  Hyades,  an  age  t  of  roughly  600  Myr  does  the  least  vio¬ 
lence  to  the  masses  and  absolute  V  magnitudes  that  we  derive. 
However,  even  isochrones  near  this  age  result  in  stars  whose 
masses  are  too  small  for  their  luminosities.  For  stars  of  2. 31  and 
2.01  M. q,  these  magnitudes  are  somewhat  brighter  than  those 
expected  from  evolutionary  models  fit  to  other  Hyades  binaries. 
The  disagreement  may  not  be  surprising  in  view  of  the  difficul¬ 
ties  in  modeling  near  the  main-sequence  turnoff. 

For  instance,  the  point  where  the  7  =  0.28,  Z  =  0.026, 
t  =  700  Myr  isochrone  of  Lebreton  et  al.  (2001)  passes  through 
the  color  and  magnitude  of  the  secondary,  B  —  V  =  0.18  mag, 
Mv  =1.61  mag,  corresponds  to  a  mass  of  1 .98  M  0,  just  within 
our  68%  confidence  contour.  But  this  isochrone  does  not  pass 
near  the  point  B  —  V  =  0.18  mag,  My  =  0.48  mag  that  corre¬ 
sponds  to  the  primary  of  62  Tau.  At  B  —  V  =  0.18  mag,  this 
isochrone  has  M  =  2.28  M  Q,  an  acceptable  mass,  but  it  is  too 
faint  at  My  =  0.8  mag.  On  the  other  hand,  if  we  boost  the  lu¬ 
minosity  to  set  My  =  0.49  mag  (the  nearest  tabulated  point),  the 
isochrone  has  the  component  as  too  massive  (A4  =  2.41  M  Q) 
and  too  red  {B  —  V  =  0.26  mag).  Older  isochrones  are  too  red  to 
match  the  secondary. 

The  7  =  0.30,  Z  =  0.030,  t  =  560  Myr  isochrone  of  Girardi 
et  al.  (2000)  fares  slightly  better,  as  illustrated  in  Figure  3.  It 
passes  through  the  color  and  magnitude  of  the  secondary  at  a 
mass  of  2.00  A4  0,  and  through  the  color  and  magnitude  of  the 
primary  at  a  mass  of  2.44  A4g.  This  pair  of  masses  is  just  be¬ 
yond  the  90%  confidence  contour.  The  isochrones  of  Lejeune  & 
Schaerer  (2001)  and  the  empirically  adjusted  Hyades  isochrone 
of  Pinsonneault  et  al.  (2004)  show  similar  results.  (In  Fig.  3a, 
the  more  heavily  marked  parts  of  the  isochrones  correspond  to 
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Fig.  2. — Confidence  contours  of  the  component  masses  (left) 
confidence  regions. 


and  absolute  V  magnitudes  (right)  of  02  Tau.  Contour  levels  are  for  the  50%,  68%,  and  90% 


the  1  a  ranges  of  allowed  masses  for  the  two  components,  and 
the  points  corresponding  to  the  best-fit  masses  are  marked  with 
filled  circles.) 

The  agreement  between  the  isochrones  and  our  results  would 
be  considerably  better  if  Z  took  the  solar  value  of  0.0 19,  although 
this  value  would  not  fit  the  rest  of  the  Hyades.  The  gray  line 
in  each  panel  of  Figure  3  is  an  isochrone  at  /  =  750  Myr,  Z  = 
0.019,  interpolated  between  the  710  and  790  Myr  isochrones  of 
Girardi  et  al.  (2000).  The  observed  9  2  Tau  B  color  and  magnitude 
fit  this  isochrone  well  along  a  stretch  of  the  isochrone  that  also 
corresponds  to  the  correct  mass  within  its  uncertainties.  This  iso¬ 
chrone  passes  through  the  92  Tau  A  color  and  magnitude  within 
their  uncertainties,  but  at  the  wrong  mass.  However,  the  error  bar 
in  Figure  3  b  shows  that  the  mass  discrepancy  for  component  A 
is  barely  more  than  1  a.  An  isochrone  of  Lejeune  &  Schaerer 
(2001)  with  t  =  790  Myr  and  solar  metallicity  is  similarly  close 
to  the  best-fit  values  of  color,  magnitude,  and  mass. 

The  situation  is  somewhat  unsatisfactory:  the  isochrones  that 
fit  are  only  a  little  older  than  current  estimates,  but  they  are  sig¬ 
nificantly  metal-poor  for  a  member  of  the  Hyades.  To  get  metal- 
rich  isochrones  that  agree  with  our  results,  one  would  have  to 
increase  the  helium  content  significantly,  which  is  at  odds  with 
the  constraints  placed  on  Tby  the  low-mass  Hyades  double-lined 
eclipsing  binary  vB  22  by  Lebreton  et  al.  (2001),  0.255  ±  0.006, 
and  Pinsonneault  et  al.  (2003),  0.271  ±  0.006,  both  less  than  the 
solar-scaled  value  of  0.28. 

6.2.  Effects  of  Rotation 

Can  the  discrepancies  be  resolved  by  invoking  the  effects  of 
stellar  rotation  on  the  observed  colors  and  magnitudes,  discussed 
in  the  context  of  Altair  by  Peterson  et  al.  (2006)?  Our  results 

TABLE  5 

Component  Masses  and  Absolute  Magnitudes 


Parameter  A  B 

Mass  (M  e) .  2.15  ±0.12  1.87  ±  0.11 

Mv  (mag) .  0.48  ±  0.05  1.61  ±  0.06 


B  -  V 


Jl  (JIQ) 

Fig.  3. — Components  of  02  Tau  compared  to  Girardi  et  al.  (2000)  iso¬ 
chrones  with  Y  =  0.30.  For  the  colored  isochrones,  Z  =  0.03.  An  isochrone 
with  Z  =  0.019,  in  gray,  is  shown  for  comparison.  ( a )  Color-magnitude  dia¬ 
gram.  The  thicker  portions  of  the  isochrones  mark  the  ±  1  a  mass  ranges.  The 
best-fit  masses  are  marked  by  filled  circles.  Crossed  black  error  bars  mark  the 
observed  B  —  V  and  My  for  each  component.  Open  circles  with  crosses  show 
the  locations  of  nonrotating  stars  of  the  same  mass  and  luminosity,  for  three 
inclination  angles.  In  the  order  shown  by  the  arrows,  i  =  90°  (equator-on), 
i  =  47?6  (=z0rbX  ar>d  i  —  13?5  (rotation  at  breakup  velocity)  for  the  primary; 
and  i  =  90°,  47?6,  and  1 7?3  for  the  secondary  (using  z;sinz  =110  km  s-1). 
( b )  Mass-magnitude  diagram.  Isochrones  are  color  coded  as  in  (a).  Filled  circles 
mark  where  B  —  V  —  0. 18  mag  (with  no  correction  for  rotation).  On  the  630- 
710  Myr  isochrones,  the  bluest  points  before  the  turnoff  (B  —  V  —  0.184,  0.21, 
and  0.22  mag,  respectively)  are  marked  with  open  circles.  Crossed  black  error 
bars  show  the  observed  masses  and  Mv.  Open  circles  with  crosses  show  the 
locations  of  the  nonrotating  equivalents,  in  the  same  order  as  in  (a). 
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TABLE  6 

Effects  of  Rotation 


Rotating  Star  Nonrotating  Equivalent 


Inclination  7p0lar  neq  reff 

(deg)  (K)  (km  s-1)  ui  Mv  B  —  V  (K) 

(1)  (2)  (3)  (4)  (5)  (6)  (7) 


Primary  (M  =  215  A40,  v  sin  i  =  65  km  s  ’) 

90 .  7957  65  0.35  0.46  0.17  7915 

47.6  (=iorb) .  7988  89  0.47  0.49  0.17  7910 

13.5 .  8748  300  0.999  1.03  0.20  7820 


Secondary  (M  =  187  A40,  v  sin  i=  110  km  s  ') 

90 .  8291  110  0.47  1.57  0.16  8205 

47.6  (=iorb) .  8396  149  0.561  1.62  0.16  8230 

17.3 .  9039  370  0.999  2.13  0.19  8080 


Notes. — Col.  (1):  Assumed  inclination  for  each  case.  Cols.  (2)  and  (3):  Derived  polar  effective  temperature 
and  rotational  velocity  for  the  observed  star  with  the  given  inclination.  Col.  (4):  Rotational  velocity  of  the  star 
as  a  fraction  of  breakup  velocity.  Cols.  (5)— (7):  Absolute  F  magnitude,  B  —  V  color,  and  effective  temperature 
for  the  equivalent  nonrotating  star  of  the  same  mass. 


show  that  rotation  does  not  alter  the  appearances  of  the  com¬ 
ponents  in  a  way  that  resolves  these  discrepancies. 

Rapid  rotation  has  several  effects.  For  stars  in  this  mass  range, 
the  luminosity  and  polar  radius  decrease  slightly  (e.g.,  Sackmann 
1970)  while  the  equatorial  radius  increases  significantly,  reach¬ 
ing  3/2  the  polar  radius  at  breakup  angular  velocity  (Hardorp  & 
Strittmatter  1968).  The  effective  temperature  of  the  equatorial 
region  decreases  significantly;  in  the  case  of  Altair,  an  A  star 
with  a  mass  of  1 .8  Ai  0  and  rotating  at  90%  of  breakup  velocity, 
the  pole-to-equator  temperature  difference  is  1850  K  (Peterson 
et  al.  2006). 

Thus,  a  rotating  star  does  not  radiate  isotropically  (Kraft 
1970).  Because  of  the  variation  in  intensity  across  the  surface,  a 
star  rotating  near  breakup  velocity  and  seen  nearly  pole-on  ap¬ 
pears  significantly  brighter  than  it  would  if  it  were  not  rotating, 
in  spite  of  the  slight  reduction  in  total  luminosity.  A  star  with  the 
same  v  sin  i  but  rotating  equator-on  has  a  smaller  rotational  ve¬ 
locity;  thus,  its  luminosity  is  only  slightly  reduced  and  the  pole- 
to-equator  variation  is  less.  The  net  effect  in  this  case  is  that  the 
star  appears  slightly  fainter  and  somewhat  redder  than  the  non¬ 
rotating  equivalent,  since  our  view  of  the  star  is  dominated  by  the 
cooler  equatorial  regions.  Figueras  &  Blasi  (1998)  have  shown 
that  the  ages  of  open  clusters  derived  from  upper-main-sequence 
isochrone  fitting  can  be  overestimated  by  30%-50%  if  the  ef¬ 
fects  of  rotation  on  ubvyHB  colors  are  not  accounted  for. 

We  used  the  suite  of  programs  developed  by  Peterson  et  al. 
(2006)  to  calculate  the  run  of  specific  intensity  across  the  surface 
of  a  Roche  spheroid  and  the  resulting  appearance  of  such  a  star 
as  seen  from  a  specific  direction.  The  models  were  computed  at 
the  masses  determined  in  this  work  and  were  constrained  to  re¬ 
produce  the  observed  V magnitude,  B  —  V  color,  and  v  sin  i.  They 
were  calculated  for  three  inclination  angles  for  each  of  the  binary 
components:  at  the  minimum  possible  inclination,  i.e.,  the  in¬ 
clination  that  implies  that  the  star  is  rotating  at  breakup  angular 
velocity;  at  i  =  i0T t,;  and  at  i  =  90°,  i.e.,  equator-on. 

Figure  3  shows  the  results  of  our  calculations  as  open  circles 
with  crosses:  the  upper  three  points  for  the  primary  at  2.15  Ai 0 
and  the  lower  three  for  the  secondary  at  1.87  A4  0.  Each  point 
shows  the  absolute  V  magnitude  and  the  B  —  V  color  of  the 
nonrotating  equivalent  star  of  the  appropriate  mass  that,  if  it  were 
rotating,  would  match  the  observed  values  of  V,  B  —  V,  and 


v  sin  i.  The  i  =  90°  and  i  =  zorb  cases  produce  nearly  identical 
values  for  My  and  B  V  for  the  nonrotating  equivalents,  while 
for  the  minimum  inclination,  indicated  by  the  arrow,  the  non¬ 
rotating  equivalent  stars  are  considerably  fainter  and  somewhat 
redder.  The  results  are  also  summarized  in  Table  6. 

The  results  of  these  calculations  do  not  substantially  improve 
the  agreement  between  our  observations  and  the  evolutionary 
models.  To  demonstrate  this  point,  we  replot  in  Figure  4  the 
models  shown  in  Figure  3,  but  as  evolutionary  tracks  for  Ai  = 
2.21,  2.15,  and  2.03  Ai 0  (the  best-fit  mass  [±o\vf])  for  the  pri¬ 
mary,  and  for  M  =  1.98,  1.87,  and  1.76  A4@  for  the  secondary. 
The  nonrotating  equivalents — calculated  at  2. 1 5  and  1 .87  Ai  0 — 
are  repeated  from  Figure  3  as  well. 

For  the  primary,  the  2.15  Ai0  nonrotating  equivalent  falls 
near  the  2.15  Ai0  evolutionary  track  near  an  age  of  670  Myr, 
but  only  if  the  primary  itself  is  rotating  at  breakup  velocity.  For 
the  secondary,  the  situation  is  worse:  the  1.87  Ai0  nonrotating 
equivalent  crosses  the  1.87  At  0  track  considerably  blueward  of 
the  500  Myr  point,  the  youngest  one  plotted,  and  then  only  if  the 
secondary  itself  is  rotating  close  to  breakup  velocity. 


B  -  V 

Fig.  4. — Nonrotating  equivalents  of  the  components  of  62  Tau  compared  to 
evolutionary  tracks  from  Girardi  et  al.  (2000)  with  Y  —  0.30.  Tracks  are  plotted 
for  stars  with  masses  equal  to  those  of  the  primary  and  secondary,  .  The 
colors  and  magnitudes  of  the  2.15  and  1.87  A't©  nonrotating  equivalents,  i.e., 
nonrotating  stars  that,  were  they  rotating,  would  produce  the  observed  V,B  —  V, 
and  v  sin  i,  are  also  plotted. 
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7.  SUMMARY 

We  have  measured  the  visual  orbit  of  the  spectroscopic  binary 
0 2  Tauri  with  the  Mark  III  optical  interferometer  and  the  Navy 
Prototype  Optical  Interferometer.  The  results,  combined  with 
other  data  in  the  literature,  allow  us  to  determine  the  masses  and 
absolute  magnitudes  of  the  components  and  to  test  some  of  the 
current  stellar  evolution  results. 

We  used  our  measurement  of  the  angular  size  of  the  orbit  in 
conjunction  with  the  parallax  of  de  Bruijne  et  al.  (2001)  and  the 
orbital  period  of  Torres  et  al.  (1997)  to  calculate  the  total  mass 
of  the  system  by  way  of  Kepler’s  third  law.  In  order  to  choose 
among  the  three  spectroscopic  orbits  in  the  literature,  we  used 
each  of  them  to  estimate  the  total  mass  and  chose  the  one  (Torres 
et  al.  1997)  that  agreed  best  with  the  parallax  result.  We  then  used 
the  ratio  of  velocity  amplitudes  to  allocate  the  total  mass  between 
the  two  components. 

The  resulting  masses  are  1.87  ±  0.11  and  2.15  ±  0.12  M& 
considered  individually;  considered  jointly,  masses  of  2.01  and 


2.31  M.q  lie  just  within  the  68%  confidence  contour.  The  abso¬ 
lute  V magnitudes  are  0.48  and  1.61  mag.  The  estimated  masses, 
especially  the  mass  of  the  primary,  are  somewhat  small  compared 
to  those  implied  by  the  observed  absolute  magnitudes  and  evo¬ 
lutionary  models  for  the  age  and  metallicity  of  the  Hyades.  We 
also  calculated  the  possible  effects  of  rotation  on  the  magnitudes 
and  colors  of  the  components.  The  results  do  not  reconcile  the 
measured  masses  with  the  evolutionary  models. 

We  thank  Lu  Rarogiewicz  for  years  of  unfaltering  assistance 
in  operating  the  Mark  111  interferometer,  Craig  Denison  for  even 
more  years  of  careful  and  conscientious  assistance  with  both  the 
Mark  Ill  and  the  NPOl,  and  H.  R.  Schmitt,  R.  T.  Zavala,  B.  D. 
Mason,  and  H.  C.  Harris  for  careful  reading  of  the  manuscript. 
This  work  was  funded  by  the  Office  of  Naval  Research  and  the 
Oceanographer  of  the  Navy.  We  have  made  use  of  the  S1MBAD 
database  operated  by  the  CDS,  Strasbourg,  France,  and  of  NASA’s 
Astrophysics  Data  System. 
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